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Abstract 

A general shell model formalism for the nonmesonic weak decay of the hypernuclei has been 
developed. It involves a partial wave expansion of the emitted nucleon waves, preserves naturally 
the antisymmetrization between the escaping particles and the residual core, and contains as a 
particular case the weak A-core coupling formalism. The hypernuclei are grouped having in view 
their A — 1 cores, that is in those with even-even, even-odd and odd-odd cores. It is shown 
that in all three cases the nuclear structure manifests itself basically through Pauli Principle, 
and very simple expressions are derived for the neutron and proton induced decays rates, T n and 
T p , which does not involve the spectroscopic factors. For the strangeness-changing weak AN — > 
NN transition potential we use the One-Meson- Exchange Model (OMEM), which comprises the 
exchange of the complete pseudoscalar and vector meson octets (n, n, K, p, uj, K*). We evaluate 3 H, 
A H, \He, \He, ll B, ^ 2 C, O, \ 7 0, and j^Si hypernuclei, with commonly used parametrization for 
the OMEM, and compare the results with the available experimental information. The calculated 
rates T^m = r n + T p are consistent with the data, but the measurements of T n / p = T n /T p are 
not well accounted for by the theory. It is suggested that, unless additional degrees of freedom are 
incorporated, the OMEM parameters should be radically modified. 



I. INTRODUCTION 



Hypernuclear physics adds another flavor (strangeness) to the traditional nuclear physics, 
and its goal is to study the behavior of hyperons (A, E, S, Q) in the nuclear environments, 
which are now bound system of neutrons, protons and one or more hyperons. Interesting 
strange nuclei with strangeness S — — 1 are the A hypernuclei, in which a A hyperon, having 
a mass of 1116 MeV and zero charge and isospin, replaces one of the nucleons. Same as the 
free A hyperon, they are mostly produced via the strong interactions, i.e., in the reaction 
processes 7i + n — > AK + , K~n — > n~A and K~p — > 7r°A, by making use of the pion (ir) and 
kaon (K) beams. They also basically decay through the weak interactions, as the free A 
does. Yet, as it is well known and explained below, there are some very important differences 
in the corresponding decaying modes. 

First, it should be remembered that the free A hyperon decays nearly 100 % of the time 
by the A — > Ntt weak-mesonic mode (Fig. 1): 




A 



FIG. 1: Mesonic (nonleptonic) decay vertex T~0^ NlT . 
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p + TT~ (64.1%) 

n + 7T° (35.7%), 

with the total transition rate + r° = T = 2.50 • 10 -6 eV (which corresponds to the 
lifetime r° = 2.63 • 10~ 10 sec). For the decay at rest the energy-momentum conservation 
implies 



MA = M N + 7 ^- + ^pl + ml; 



PN = Ptt- 



Therefore the energy released is 



Q M A M N ///, ~~ 37 MeV, 

and the kinetic energies and momenta in the final state are: 

T N = (MA ~ ^ )2 + m ^ 5MeV; T n = Q -T N = 32 MeV, 
ZM A 

Pn = Pn = y/{T N + M N ) 2 - M 2 N = 100 MeV/c. 

It is clear that the isospin is changed by AT = 1/2 and 3/2 and its projection by AM T 
— 1/2 in the free A — > Nix decay. More, as 

• A -»■ n + 7T°: 

\n)\n°) = |l/2, -1/2)|10) = 7||l/2, 1; 1/2, -1/2) + y||l/2, 1; 3/2, -1/2) 



A — > p + 7r : 

|p>|0 = |l/2, 1/2, )|1, -1) = - j||l/2, 1; 1/2, -1/2) + j||l/2, 1; 3/2, -1/2), 



one sees that the above experimental data can be accounted for fairly well by neglecting 
the 1 1/2, 1; 3/2, — 1/2) components in these relations (AT = 1/2 rule). In fact, one gets 
T^-fT^o = 2, while the experimental result is 64.1/35.7 = 1.80. 

Assuming the AT = 1/2 rule, the phenomelogical weak Hamiltonian for the process 
depicted in Fig. 1 can be expressed as: 

HZfn = -^G F m 2 J N (A n + B nl5 ) <j> w ■ t^a ^ j , (1.1) 

where Gpfn 2 = 2.21 x 10~ 7 is the weak coupling constant. The empirical constants A n = 1.05 
and B n = —7.15, adjusted to the observables of the free A decay, determine the strengths 
of parity violating and parity conserving amplitudes, respectively. The nucleon, A and 

pion fields are given by ipN and ipA and <f> w , respectively, while the isospin spurion 

1 



is included in order to enforce the empirical AT = 1/2 rule. (Note that: iPn4>-k ' — 

(V>jAh> + V2tp n <j)n- , V / 2^p07r+ ~ A). 

The free A hyperon weak decay is radically modified in the nuclear environment because 
the nucleon and the hyperon now move, respectively, in the mean fields Un and U\, which 
come from the NN and iVA interactions. Un and U\ are characterized by the single particle 
energies (s.p.e.) En and e\ and we have to differentiate between: 

1. Mesonic Decay (MD): The basic process is again represented by the graph shown in 



Fig. 1 and described by the hamiltonian ( pL.ip . Yet, the energy-momentum conservation is 
different: 



M A = M N -s A + s N + + \jp\ + ml ; p A = - PlT 

where Ma = AMn and are, respectively, the mass and the momentum of the whole 
nucleus; A is the mass number, and e N are the s.p.e. of the loosely bound states above the 
Fermi energy e N . They are of the order of a few MeV, while E\ is the energy of the Osi/2 
state and goes from —11.7 MeV for A 3 C to —26.5 MeV for A 08 Pfe Thus, the corresponding 
Q- values 



Qm .U a M x - //>, A 



T 

N " 



are significantly smaller than Q , particularly for medium and heavy nuclei. The experi- 
mental decay rates T n - + r^o = Tm = Tm(A — > Nn) are of the order of T° only for nuclei 
with A < 4, and they rapidly fall as a function of nuclear mass. For instance, in A 2 C: 
r^o/r = 0.217 ± 0.084 and r^-/r° = 0.052l£$|. This hindrance effect, as is illustrated in 
Fig. 2 for the hypernucleus A 7 0, is due to Pauli principle. In fact, the population of states 
that are below the Fermi level, with energies e N < e N , is totally blocked by the Pauli prin- 
ciple, while transitions to states that lie above e N are strongly hindered due to the selection 
rule AN= 0, N being the harmonic oscillator quantum numbers. That is, only a few second 
forbidden transitions (AN= 2) can occur in the case of A 7 0, and it is clear that the degree 
of forbiddiness increases with A. 

2. Nonmesonic Decay: New nonmesonic decay (NMD) channels AA^ — > NN become 
open inside the nucleus, where there are no pions in the final state. The corresponding 
transition rates can be stimulated either by protons, T p = T(Ap — > np), or by neutrons, 




FIG. 2: Schematic picture of the hypernuclear mesonic decay in \ 7 0: The population of occupied 
states 0s 1 / 2 ,0p 3 / 2 and 0p 1 / 2 , which are below the Fermi level £ N , is totally blocked by the Pauli 
principle, while transitions to the weakly bound empty states 0d$/2 and ls^, which lie above e N , 
are strongly hindered by the selection rule AN= 0. 



T n = T(An — > nn). The energy- momentum conservation and the Q- value are, respectively: 

2 2 2 

M A = M N - e A - e l N + ^- + + p, = p - p .. 
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and 



Qnm = M A - M N + e A + e l N , 



where pi and p 2 are the momenta of the two outgoing nucleons. As the mean value of ejy is 
~ 30 MeV one gets that Qnm ~ 120 — 135 MeV, which is basically the kinetic energy of the 
two particles that are ejected from the hypernucleus. This means that the nonnesonic decay 
process possesses a large phase space in the continuum, as is outlined in Fig. 2 for the case of 
the hypernucleus A 0. The theoretical models reproduce fairly well the experimental values 
of the total width T NM = T n + T p (F e ^ M = F°) but the ratio T n/p = T n /T p (0.5 < r^P < 2) 
remains a puzzle. 




FIG. 3: Schematic picture of the hypernuclear nonmesonic decay in j^O: The hyperon A and one 
of the nucleons from the occupied states 0si/ 2 ^P3/2 an d Opi/2 are expelled into the large phase 
space in the continuum, becoming free nucleons with momenta pi and p2- 



Very often it is assumed that the hypernuclear NMD AiV — > A^ A" is triggered via the 
exchange of a virtual meson, and the obvious candidate is the one-pion-exchange (OPE) 
mechanism, where the strong Hamiltonian 



n NN T , 



igNN-K^NlsK ■ Tlf) N , 



[1.2) 



(with Qnn-k = 13.4) accompanies the weak Hamiltonian (LI). Following the pioneering 
investigations of Adams [13| several calculations have been done within this coupling scheme 
yielding: Y^JP = T° and T^ E) = 0.1 - 0.2 @, [T5J, [16], |I7|, 0, |Tj| |20], [21], 0, pj, ^ |25|, 
ESI, E7L BR E9l BUI RTl B3, B3L B|, RR BE B7L B3, 531 . The importance of the p meson in the 



weak decav mechanism was first discussed bv McKellar and Gibson HT31, and the present- 
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FIG. 4: The nonmesonic weak decay mode A + N — > N + N. In the one meson exchange model 
a weak vertex Tiw is always combined with a strong vertex Tis- 

day consensus is, however, that the effect of the p-meson on both T NM and T n / p is small 
Jj| |27|, |3£|. The full one meson-exchange model (OMEM), which encompasses all 



pseudoscalar mesons (ir, rj, K) and all vector mesons (p, u, K*), has been also considered by 
several authors |]24], |5|, |34], [35], [36[ |3^, ^9[. From these works we have learned that, although 
the K meson contribution significantly increases the ratio T n / P , the OMEM is unable to 
account for the corresponding experimental values. 

We wish to restate that the OMEM transition potential is purely phenomenological and 
that it is not derived from a fundamental underlying form, as happens for instance, in the 
case of electro-magnetic transitions or the semileptonic weak decays. As in the case of the 
OPE, in the OMEM, a weak baryon-baryon-meson (BBM) coupling is always combined 
with a strong BBM coupling (see Fig. 4). The strong one is determined experimentally with 
some help from the SU(3) symmetry, and the involving uncertainties have been copiously 
discussed in the literature j|D], ^T], f£|. It is the weak BBM couplings which could become 
the largest source of errors. In fact, only the weak AA7r amplitude can be taken from the 
experiment, at the expense of neglecting the off-mass-shell corrections. All other weak BBM 
couplings are derived theoretically by using SU(3) and SU(6) W symmetries, octet dominance, 
current algebra, PCAC, pole dominance, etc. pi M, M, g|, [45], g6], |7[ ||, M, M, [51 



Assortments of such methods have been developed and employed for a long time in weak 
interaction physics to explain the hyperon nonleptonic decays. One should also keep in 
mind that both the strong and weak BBM couplings, as well the meson masses, can become 



significantly renormalized by the nuclear environment [52] 



The nuclear structure frameworks utilized in the literature for the formal derivations of 
the NMD rates are: i) the nuclear matter, and ii) the nuclear shell model. There are relatively 
few works where the second method was employed, and up to quite recently all they were 
involved the technique of coefficients of fractional parentage, with the spectroscopic factors 



(SF) explicitly appearing in the expressions for the transition rates |20], g^, |30], [38]. At 
variance, Barbero et al. |39j have developed a fully general shell model formalism and have 
specified it for hypernuclei with odd-mass core, such as \H, \He, jfC, and j^Si, were the 
cores are: 3 H, 3 He, U C and 21 Si. Here we also discuss the hypernuclei with even-even and 
odd-odd number of protons and neutrons, namely \He and j^O, and \H, and X ^B. 



II. SHELL MODEL FORMALISM 



The shell model framework for the NMD rates has been developed in detail in Ref. |3l 
and here we will just sketch the main steps, employing the same notation. One starts from 
Fermi's Golden Rule for the decay rate, 

r = 2tt Y, I |(PiP 2 W 5 ,J F M F ;TM T |y|J / M 7 )| 2 

SM s JfM f TM t j 

dpi dp 2 



x 5(e Pl +e P2 +£ F -S! 



(2.1) 



(27r) 3 (27r) 3 ' 

where V is the weak hypernuclear potential, and the wave functions for the kets 
\p1P2SMs, JfMf; TMt) and \JjMi) are assumed to be antisymmetrized and normalized. 
After performing: 1) the transformation to the relative and center of mass (cm.) momenta, 
p and P, and angular momenta 1 and L, and 2) the angular momentum couplings: 1 + L = A, 
A + S = J one obtains: 



16M 3 



7T 



N JJ 2 



E 



SXlLTJJ F a 



deJe(A% - e) 



£ M(pPlLXSJT;j A j N M)(Ji\\ (°<] N t N alt A ) jK) 



3Nlh 



(2.2) 



where J = y/2J + l, P = 2^/M^, p = ^ M N (A - e), A£ = E 1 - £f - 2M N , t A = -1/2, 
t p = 1/2, t n = —1/2, and the label a goes over all final states with the same spin and 



parity. Thus, that the NMD rates, in principle, depend on both: i) the nuclear structure 
effects through the two-particle NA parentage coefficients (Ji\ \ (&] N t N a j A tA) , W^f), an d on 
the transition potential via the elementary transition amplitudes 

M(pPlLXSJT;j AjN ,t N ) - ' 1 ' ''- s - v 'i 



> 



x (pPlL\SJ-TM T = t A + t N \V\ 3A3N J;t A t N ). (2.3) 

Here (■ ■ • |V| ■ ■ •) is the direct matrix element and the factor in front takes care of the 
antisymmetrization. To evaluate the nuclear matrix element one has to carry out the jj — LS 



recoupling and the Moshinsky transformation |53| on the ket \j A jNJ)'- 



\jaJnJ) = JaJn E 

A'5'nlNL 



In \ 3n 
X' S' J 



> (n\Nl\'\n A l A n N l N X')\n\Nl\'S'J), (2.4) 



where (•••!•••) are the Moshinsky brackets ]5B , and I and L stand for the quantum numbers 
of the relative and cm. orbital angular momenta in the AN system. The explicit expressions 
for the transition potentials are given in the Ref. |39|| . 

When the hyperon is assumed to be weakly coupled to the A — 1 core, which implies that 
the interaction of A with core nucleons is disregarded, one has that \Jj) = \(Jc3a)Ji), where 
Jc is the spin of the core, and gets 



(Ji\\(a) 



3N^N 3A^A J J 



\Jf) 



Jc Ji 3a 
J j N J F 



,j f+ j +Ji jj i{ ^ - - \(Jc\\a) Nt JJ F ). (2.5) 



To evaluate the one-particle spectroscopic amplitudes {Jc\Wj N t N \\Jp) we will use the BCS 
approximation and, it will be assumed that the even-even, odd-even and odd-odd cores are 
described, respectively, as zero, one and two quasiparticle states. Correspondingly, the state 
\J C ) goes into \BCS), b) it \BCS) and (b] lh b] 2t2 )j 2 \BCS}, where \BCS) is the BCS vacuum 



and bj = UjOj — Vjaj is the quasiparticle creation operator In all three cases the NMD 
rate can be cast in the form: 



J=3N+3A 

£ 

On J=\3n~3a\ 



T t N = £ £ tf N Fj(j N t N )Kj(j N t N ) } 



where 



N, 



7T JO 



SILXT 



and 



Aj N t N - M A - M N + e jAtA + e jNtN . 
The geometrical factors Fj(j N t N ) are: 

1. Even-Even core: \Jc) — ► \BCS), 



(2.6) 



v / de^e(A jNtN - e) £ M 2 (pPlLXSJT; j A j N ,t N ), (2.7) 



(2.8) 



Fj{j N t 



A' J 



3a 2 J 2 



(2.9) 



2. Odd-Even core: \J C ) -»• b] ltl \BCS), 



31+0 N 

Fj(j N t N ) = J 2 J2 

Jf=\3i-0n\ 



I + (—) Jf 8 jljN 5 tNtl 



92 

•J T? 



ji 3n Jf 
J Ji 3k 



(2.10) 



3. Odd-Odd core: \J n ) -> (&L. &L_ ) d BCS) , 



Jl+32 jl+3N 

Fj(j N t N ) = 5 tNtl J 2 J 2 2 Jl £ Jl [l + (-) Jl 5 jljN 8 tNtl 

J F =\Jl-j z \ Jl=\jl—3N\ 
2 

+ (jlh) <"> C/2^) 
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' Jf 


Jl 


■M 
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JN 






J A 


JN 





It is clear that in (|2.11 ) £ 2 = — t 



(2.11) 



TABLE I: Values of (2jjv + l)Fj(jj^t]y). The quantum numbers (ji,ii, J/) in odd-even core nu- 
clei are: (0s 1/2 , -1/2, 0) for ^Te, (0s 1/a , 1/2,0) for \H, (0p 3/2 , -1/2, 1) for \ 2 C, (0p 1/2 , -1/2, 1) 
for A e O, and (0cZ 5 /2, — 1/2, 2) for ^S 1 *. Similarly, (ji, t\, ji, t\; J2J1) in odd-odd core nuclei are: 
(0s 1/2 , -1/2, 0s 1/2 , 1/2; 1,1/2) for \He, and (0p 3/2 , -1/2, 0p 3/2 , 1/2; 3, 5/2) for ^B. 
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The nuclear structure manifests itself basically through the factors Fr(iivtiv), which are 



engendered by the Pauli principle. Their values for a few cases are given in Table I. We 
would like to stress that the quantum numbers jit\ and stand for the hyperon partners 
in the initial state, and that runs over all proton and neutron occupied states in the 
initial nucleus. It is amazing to notice that the the last three equations are valid for any 
hypernucleus, which could be so light as \H or so heavy as f 8 Pb. One should also add 
that the Eq. ( [2.6|) contains the same physics as the Eq. (5) in Ref. [25[] or the Eq. (30) in 



Ref. ||38|| , with the advantage that we do not have to deal with spectroscopic factors. Of 
course, neither the initial and final wave functions are needed. From the results displayed in 
Table I it can be seen that in all three cases the coefficients Fj(j '^t n) are of the same order 
of magnitude, which indicates that the nuclear structure effects in the NMD are of minor 
importance. 



III. NUMERICAL RESULTS AND DISCUSSION 



The numerical values of the parameters necessary to specify the transition potential, were 
taken from Ref . E5J , where, in turn, the strong couplings have been taken from Refs. [HO, 



and the weak ones from Ref. [23|. The energy difference Aj N t N i n (^) is evaluated from 
the experimental single nucleon and hyperon energies. It is a general belief nowadays that, 
in any realistic evaluation of the hypernuclear NMD rates, the finite nucleon size (FNS) and 
the short range correlations (SRC) have to be included simultaneously. Therefore, in the 
present paper both the FNS and SRC renormalization effects are considered, in the way 



described in Ref. |39| . Under these circumstances, and because of the relative smallness of 



pion mass, the transition is dominated by the OPE [25] , [3£ 



The numerical calculations were done in the extreme shell model, which implies that the 
pairing factors Vj N were taken to be equal to one (zero) for the occupied (empty) levels. 
Thus, from the nuclear structure point of view, the only free parameter is the harmonic 



oscillator length b. We evaluate it from the relation 6 = 1/ \/KujMn, and the oscillator 
energy was estimated from the relation %u = — 25A~ 2 / 3 MeV, which is frequently 

used for light nuclei. The calculations for T^m and T n / p are confronted the experimental 
data in Tables [TT] and in |TJ, respectively. Three different OMEM have been employed for the 
transition potential. Namely: (tc) only the pion was taken into the account, (PS) all three 
pseudoscalar mesons (tt + r\ + K) were included, and (PS + V) also the vector (p + uj + K*) 
mesons are considered. The same remarks are pertinent here as in the study |39| where only 
j^C has been analyzed. That is: (1) the simple OPE model accounts for Tnm, but it fails 
badly regarding T n / P , (2) when r\ and K mesons are included, the total transition rate is only 
slightly modified, while T n / p change significantly, coming somewhat closer to the measured 
values, and (3) the results are not drastically modified when all vector mesons are built-in. 



IV. SUMMARY AND CONCLUSIONS 



The shell model formalism for the nonmesonic weak decay of the hypernuclei involves a 
partial wave expansion of the emitted nucleon waves and preserves naturally the antisym- 
metrization between the escaping particles and the residual core. The general expression 
( p.2|) is valid for any nuclear model and it shows that the nonmesonic transition rates should 



TABLE II: Nonmesonic decay rates Tnm m units of T° = 2.50 • 10 eV. The symbols tt, PS 
and V stand, respectively, for the transition potentials activated by the pion, the pseudo-scalar 
(tt + r] + K), and the vector (p + u + K*) mesons. Experimental data are shown for comparison. 
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fSi 


1.314 


1.100 
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1.30 ±0.10 0] 



depend, in principle, on both: (i) the weak transition potential, through the elementary 
transition amplitudes Ji4(pPlL\SJT; JaJn^n), and (ii) the nuclear structure, through the 
two-particle iVA parentage coefficients (Jj\\ (o] N i N a j A ) J \\Jf)- The latter explicitly depend 
on the initial and final wave functions. Yet, as explained in Ref. |ft9|| , and because of: a) the 
inclusive nature of the nonmesonic decay, and b) the peculiar properties of the coefficients 
FjUn^n)-, this dependence is washed out. In this way we have arrived at a very simple re- 
sult for transition rates, given by the Eq. ( |2.6| ), which is valid, not only for the hypernuclei 
with odd-even core (as shown previously |39]]), but also for those which have even-even and 
odd-odd cores, which has been demonstrated here. 

We reproduce satisfactorily the data for the total transition rates with the OMEM 



parametrization from the literature p5j , but the n/p-ratios are not well accounted for. 
Thus, after having acquired full control of the nuclear structure involved in the process, 
and after having convinced ourselves that the nuclear structure correlations can not play a 



TABLE III: The p/n ratios for several hypernuclei. See the Table II caption. 
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1.38 
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crucial role, we firmly believe that the currently used OMEM should be radically changed. 
Either its parametrization has to be modified or additional degrees of freedom have to be 
incorporated, such as the factorizable terms |j54| , the axial- vector- meson exchanges |55 



or 



the correlated the correlated 2n from Ref. 38 
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